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We reformulate the theory of polycrystalline plasticity, in externally driven, nonequilibrium situa- 
tions, by writing equations of motion for the flow of energy and entropy associated with dislocations. 
Within this general framework, and using a minimal model of thermally assisted depinning with 
essentially only one adjustable parameter, we find that our theory fits the strain-hardening data 
for Cu over a wide range of temperatures and six decades of strain rate. We predict the transition 
between stage II and stage III hardening, including the observation that this transition occurs at 
smaller strains for higher temperatures. We also explain why strain-rate hardening is very weak up 
to large rates; and, with just one additional number, we accurately predict the crossover to power- 
law rate hardening in the strong-shock regime. Our analysis differs in several important respects 
from conventional dislocation-mediated continuum theories. We provide some historical background 
and discuss our rationale for these differences. 
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I. INTRODUCTION 
A. Historical Background 

The time seems ripe for a critical reexamination of the 
theory of dislocation-mediated plasticity. The basic ele- 
ments of dislocation theory were described decades ago 
in the classic books by Cottrell [1| , Friedel @ , and Hirth 
and Lothe Q . Since then, major progress has been made 
via numerical simulations and experimental observations; 
but several fundamental issues remain unresolved. Most 
notably, a first-principles theory of strain hardening has 
yet to be developed. 

The question of whether such a theory is feasible re- 
mains a matter of serious debate. The prevailing opinion 
among experts in the field is that it is not. For exam- 
ple, Cottrell Q has argued recently that strain hardening 
[rather than turbulence] , is "the most difficult remaining 
problem in classical physics." He goes on to explain that 
"neither of the two main strategies of theoretical many- 
body physics - the statistical mechanical approach; and 
the reduction of the many-body problem to that of the 
behaviour of a single element of the assembly - is avail- 
able to work [strain] hardening. The first fails because 
the behaviour of the whole system is governed by that 
of weakest links, not the average, and is thermodynam- 
ically irreversible. The second fails because dislocations 
are flexible lines, interlinked and entangled, so that the 
entire system behaves more like a single object of extreme 
structural complexity and deformability, that Nabarro 
and I once compared to a bird's nest, rather than as a 
set of separate small and simpler elementary bodies." 

Similarly, on page 235 of their definitive review of 
strain hardening, Kocks and Mecking Q tell us pes- 
simistically that "an ab initio theory of strain hardening, 



with a quantitative prediction of the numerical constants, 
is unlikely to ever be derived even for a specific case, and 
impossible with any generality." They advocate what, 
from a physicist's point of view, is a purely phenomeno- 
logical approach, based on extensive observations and a 
search for trends, but with no hope of uncovering funda- 
mental principles that might lead to predictive theories. 

Devincre, Hoc and Kubin [6j are equally pessimistic 
when they say that "The present dislocation-based mod- 
els for strain hardening still have difficulties integrating 
elementary dislocation properties into a continuum de- 
scription of bulk crystals or polycrystals. As a conse- 
quence, current approaches cannot avoid making use of 
extensive parameter fitting." These authors then show 
that they can reproduce observed deformation curves by 
defining a dislocation mean free path and incorporating 
it, with input from numerical dislocation dynamics, into 
a multiscale analysis. Unfortunately, their work is based 
on a well known equation of motion for the dislocation 
density that that we find problematic on basic physical 
grounds. (See SeclIVl) 

If a first-principles theory of dislocation-mediated plas- 
ticity is intrinsically impossible, then this central part of 
solid mechanics is an intrinsically unfortunate research 
area - much less fortunate, for example, than semicon- 
ductor electronics or polymer science, where basic un- 
derstanding has been translated into quantitative theo- 
ries which, in turn, have guided enormously successful 
engineering applications. All of the authors cited above 
are keenly aware of this problem, but insist that nothing 
can be done about it. In the present paper, we challenge 
their pessimism by proposing a different approach based 
on thermodynamic principles. We examine Cottrell's ra- 
tionales for failure - especially the thermodynamically 
irreversible nature of dislocation dynamics - from this 
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broader point of view, and look to see where specific dis- 
location mechanisms play roles within a bigger picture. 



B. Thermodynamics and Dislocations 

Cottrell's remark about thermodynamic irreversibility 
reflects a general consensus that thermodynamics, espe- 
cially the second law, is irrelevant in dislocation theory, 
ft is well known that the energies of dislocations are so 
large that ordinary thermal fluctuations are completely 
ineffective in creating or annihilating them. It is also well 
known that the entropy of dislocations is extremely small 
in comparison with the total entropy of the material that 
contains them. (For example, see [lj, p. 39.) This hap- 
pens because the dislocations involve only a very small 
fraction of all the atoms in the material and, therefore, 
account for only a very small fraction of the degrees of 
freedom of the system as a whole. 

The central assertion of this paper is that, although the 
dislocation entropy is small, it is an essential ingredient 
of a theory of dislocation-mediated deformation. In writ- 
ing equations of motion for a system that contains many 
millions of irregularly moving dislocations per square cen- 
timeter, it is necessary to require that these equations al- 
ways take the system toward states of higher, and never 
lower, probability. Then, in circumstances where ordi- 
nary thermal fluctuations are irrelevant, the dislocation 
entropy by itself must be a non-decreasing function of 
time; and there must exist some dislocation-specific form 
of the second law of thermodynamics to describe this sit- 
uation. 

The natural way in which to develop a second-law anal- 
ysis for dislocations is to express the theory in terms of an 
"effective temperature," denoted here by x to distinguish 
it from the temperature T that characterizes the ordinary 
thermal fluctuations. Given a macroscopic sample of a 
material, we can compute the energy, say Uc, of any con- 
figuration of dislocations, and we can count the number 
of such configurations in any energy interval. Thus we 
can compute a configurational entropy Sc(Uc)- Follow- 
ing Gibbs, we define the effective temperature, in energy 
units, to be 



Because the dislocation energies are large, x is an ex ~ 
tremely large temperature, vastly greater than T. There- 
fore, although Sc may be small, the product X&C is of 
the order of the dislocation energy Uc', and the state that 
minimizes the free energy Uc — X&C is the most prob- 
able state of the system. An externally driven system, 
such as a deforming crystal, must be either in that most 
probable state or moving toward it. 

We emphasize that we are talking only about exter- 
nally driven systems. In the usual Gibbsian statisti- 
cal thermodynamics, one considers equilibrium states in 



which thermal fluctuations cause systems to explore sta- 
tistically significant fractions of their state spaces. If the 
dislocations are to explore their state spaces in the ab- 
sence of such fluctuations, they must be driven to do so 
by external forces. Accordingly, our thermodynamic the- 
ory of dislocations pertains only to situations in which 
there are external forces that persistently drive the sys- 
tem from one state to another, and produce complex, 
chaotic motions. It is only in such cases that it makes 
sense to enforce the second law of effective-temperature 
thermodynamics. The resulting theory is well suited for 
experiments at very high shear rates, and also should 
work for most constant strain rate experiments under lab- 
oratory conditions. But it is not relevant to quasistatic 
experiments such as indentation or beam bending; nor 
does it have anything to say (so far as we know) about 
strain-gradient theories, geometrically necessary disloca- 
tions, and the like. 

By using the word "entropy" in the same paragraph as 
the term "driven system," we raise yet another controver- 
sial issue. The question is whether the term "nonequilib- 
rium thermodynamics" makes any sense at all; or, more 
specifically, whether it is conceptually possible to derive 
the classic Kelvin-Planck or Clausius statements of the 
second law directly from Gibbsian principles of statisti- 
cal physics. A large body of literature, starting with the 
papers of Coleman, Noll, and Gurtin in the I960's 0,11], 
neatly avoids this issue by adopting the Clausius-Duhem 
inequality as an axiomatic basis for a systematic formu- 
lation of noncquilibrium thermodynamics. 

Two of the present authors (EB and JSL in 0-Qj]]) 
have found that the axiomatic approach is not sufficient 
for noncquilibrium problems in solid mechanics, where 
one must deal with internal variables such as densities 
of flow defects or dislocations, and where the internal 
degrees of freedom may fall out of equilibrium with the 
thermal reservoir. By carefully defining the entropy of an 
externally driven system containing internal variables, we 
are able to constrain the equations of motion for those 
variables to be consistent with the fundamental, statisti- 
cal statement of the second law. Our procedure underlies 
all of the following analysis. 

Finally, we emphasize that our thermodynamic formu- 
lation is not in any sense a replacement for the disloca- 
tion theories described in [lH3| and in many texts and 
research reports published since those classic works ap- 
peared. As will be seen in what follows, our analysis 
leads to the definition of a relatively small number of pa- 
rameters describing, for example, dislocation production 
rates, or the rates at which external work is converted 
to configurational entropy in various circumstances. To 
compute such quantities from first principles, we even- 
tually will need to invoke all of the specific mechanisms 
that we know to be relevant - the dynamics of different 
kinds of dislocations, their interactions with each other 
and with other internal structures, and the mechanisms 
by which they are created and destroyed. We expect 
that, in future research, we will be better able to under- 
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stand these mechanisms once they are seen in a broader 
thermodynamic context. 



C. Scheme of This Paper 

In order to make progress beyond thermodynamics, 
we need a physics-based model that relates the plastic 
strain rate and dislocation motion to the driving stress. 
In Secim we propose a minimal model of dislocation pin- 
ning and depinning that begins to address Cottrell's con- 
cerns about the irregular, irreversible, collective behavior 
of dislocations. Our model describes thermally activated 
deformation during stage-III hardening and dynamic re- 
covery, and it also provides a starting point for describing 
stage-II hardening. (See Q for a definition of the several 
"stages" of strain hardening.) However, it does not de- 
scribe the deterministically chaotic motion of dislocations 
seen at very low temperatures and extremely small strain 
rates. That behavior is discussed in detail by Zaiser [l2[ 
who, in the concluding part of his review, makes it clear 
that our range of interest does not (yet) overlap with his. 

Section [TTT] is devoted to the details of our thermody- 
namic analysis. The ideas discussed there have emerged 
from recent developments in the shear-transformation- 
zone (STZ) theory of amorphous plasticity [H El, El, 
where the key ingredient is an effective temperature that 
characterizes the internal state of disorder. The effective 
temperature in the STZ theory is exactly the same as 
the temperature \ defined in Eq. (|l.l[) . Both STZs and 
dislocations are configurational defects; that is, they are 
irregularities in the underlying atomic structures of solid 
materials. In comparison to the ephemeral STZs, how- 
ever, dislocations are long-lived, spatially extended ob- 
jects. Because dislocation energies are so much greater 
than thermal energies, and because the dislocations do 
not participate in glasslike jamming transitions and thus 
need no internal degrees of freedom of their own, the ef- 
fective temperature analysis is conceptually simpler for 
dislocation-mediated plasticity in crystalline solids than 
it is for deformations of glassy materials. We have taken 
advantage of this simplicity by writing Sec lIIII in a lan- 
guage that we intend to be self explanatory. 

A major impetus for the present investigation has been 
the paper by Preston, Tonks, and Wallace (PTW) 
on plasticity at extreme loading conditions. These au- 
thors construct a phenomenological formula for the stress 
as a function of strain, strain rate, and temperature, 
for a number of elemental metallic solids. Their for- 
mula fits direct experimental measurements and data de- 
duced indirectly from shock tests, for strain rates ranging 
from moderate to explosively large, and and for temper- 
atures ranging from cryogenic up to substantial fractions 
of melting points. It describes plastic flow as a ther- 
mally activated process for strain rates e up to at least 
10 4 sec -1 . In the explosively fast regime, for strain rates 
above about 10 8 sec -1 , the PTW curves of stress a ver- 
sus strain rate cross over to a power law, a ~ eP , with 



f3 = 0.25. The scope of the PTW analysis, plus the ap- 
parently universal nature of the phenomena that they 
describe, led us to suspect that there might be a more 
fundamental way to understand these behaviors. We use 
the PTW data for Cu in Sees I VII and I VIII for comparisons 
between our theory and experiment. 

Sections IIVI fVl and IVII contain our analysis of strain 
hardening. We start in Sec lIVI by explaining why the 
storage-recovery equation of Kocks and coworkers d, fl6l — 
fl8j is unsuitable for our purposes, and then go on to de- 
rive an alternative equation based on the second law of 
thermodynamics and energy conservation. The result- 
ing theory allows us only one fully adjustable parameter, 
specifically, the temperature (but not strain-rate) depen- 
dent fraction of the external work that is converted to 
configurational heat, as opposed to being stored in the 
form of dislocations or other internal structures. Numer- 
ical results and comparisons with experiment are pre- 
sented in Sec lVIl We show there that our theory fits the 
strain-hardening data for Cu over a wide range of tem- 
peratures and six decades of strain rate. It predicts the 
transition between stage II and stage III hardening, and 
explains why this transition occurs earlier at higher tem- 
peratures. It also explains why strain-rate hardening is 
very weak up to large rates. Finally, in Sec lVIII we show 
how a simple extension of these thermodynamic ideas 
accounts for the observed crossover to strong, power- 
law, rate hardening at extremely high strain rates in the 
strong-shock regime. [l5| 

We conclude in Sec lVIIIl with a brief summary of our 
results and remarks about future directions for investiga- 
tion. We emphasize that equations of the form proposed 
here are well suited for the study of position dependent 
instabilities such as shear banding, and we look ahead 
to comparing those equations to the ones used by Anan- 
thakrishna [19j in his studies of dynamic dislocation pat- 
terns. 



II. MINIMAL MODEL OF DISLOCATION 
DYNAMICS 

The minimal model to be used here is a polycrystalline 
material subject to a simple shear stress of magnitude 
a. Our key assumption is that dislocation motion is con- 
trolled entirely by a thermally activated depinning mech- 
anism. Thus we start with a model that we expect to be 
most relevant to stage-III hardening and dynamic recov- 
ery. We will see in Sec lIVI how this model can be used in 
a theory that includes athermal, stage II behavior. 

We assume that, after coarse-grained averaging, our 
system is spatially homogeneous and orientationally sym- 
metric; orientational symmetry is broken only by the di- 
rection of the applied stress. We further assume that we 
can describe the population of dislocations by a single, 
averaged, areal density p. For more detailed applications, 
p will need to be replaced by a set of such densities, each 
describing different kinds of dislocations with different 
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orientations, and we will need separate equations of mo- 
tion for each of these populations. Even in that case, 
however, we can assume that these populations retain 
their identities and do not, like STZs, transform dynam- 
ically from one internal state to another, thus requir- 
ing additional internal degrees of freedom. As a result, 
generalizing to a many-population version of the present 
theory should not raise fundamentally new issues. 

Suppose that the population of dislocations moves with 
average velocity v in response to a. The plastic shear rate 
is given by the Orowan relation 
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pbv; p = £ 



(2.1) 



where b is the magnitude of the Burgers vector and £ is 
the average distance between dislocations. In the spirit of 
leaving no assumption unquestioned, note that Eq. ([2.1[) 
is a purely geometric relation. Each dislocation produces 
a shear b/L as it moves across a system of linear size L. 
There are pL 2 dislocations doing this at any time; and 
the rate at which these events are occurring is v/L. The 
product of these factors is the right-hand side of Eq. (|2.1I) . 
Strictly speaking, both sides of this equation are tensors; 
but we suppress the tensor notation in interests of sim- 
plicity. 

The next step is to compute the velocity v. To do this, 
assume that each dislocation moves through an array of 
pinning sites, which may be a forest of cross dislocations. 
Alternatively, the pinning sites might be grain bound- 
aries, point defects, stacking faults, or some combination 
of all of these. Because a dislocation is a long-lived entity, 
it becomes pinned and unpinned many times during its 
motion across the system. Here we make a minor depar- 
ture from conventional dislocation theory. Rather than 
distinguishing between mobile and pinned (or "stored" ) 
dislocations, we prefer to assume that each dislocation 
retains its identity throughout these processes. 

If a dislocation spends a characteristic time Tp at a 
pinning site, and then (as in |6|) moves almost instan- 
taneously across some mean free path I* before being 
trapped again, its average speed is 

»?W = -Ct- (2-2) 
Tp{a) 

To estimate vp(a), assume a thermally activated mech- 
anism in which the dislocation is trapped in a potential 
well of depth Up(0) = ksTp in the absence of an ex- 
ternal stress, and can escape from this trap via a ther- 
mal fluctuation. The activation temperature Tp may be 
large, in fact, greater than the melting temperature; but 
it remains much smaller than the dislocation formation 
energy. When a stress a is applied, the barrier oppos- 
ing escape from the trap is lowered in the direction of a 
and raised in the opposite direction. Since we need to 
evaluate vp for arbitrarily large cr, we cannot make the 
usual linear approximation in a for this barrier-lowering 
effect, but instead must let the barrier decrease smoothly 
toward zero. The simplest way to do this is to write 



where ot is a characteristic depinning stress. 

It seems clear that ctt must be the Taylor stress: 



(2.4) 



where pr is proportional to the shear modulus p. The 
point here is that, except in the case of an isolated dis- 
location in an exceptionally pure crystal, the array of 
elastically interacting dislocations impedes the internal 
motions. The right-hand side of Eq. (l2.4p is an estimate 
of the shear stress needed to depin one dislocation seg- 
ment by moving it a distance, say b' , from its pinning 
site, where pr b/£ = pb' /£. We guess that the ratio b' /b 
is of the order of 0.1, so that pr ~ 0.1 p. On dimensional 
grounds, it is hard to see how any other combination of 
parameters could play a role in Eq. (|2.3p . 
The velocity vp(a) is 



vp (a) = — f P {<r) - /p(-ct) 
t L 



(2.5) 



where r 1 is a microscopic attempt frequency, of the or- 
der of 10 12 scc _1 , and 



fp{a) = cxp 



T 



(2.6) 



Antisymmetry is required in Eq. (|2.5[) both to preserve re- 
flection symmetry, and to satisfy the second-law require- 
ment that the energy dissipation rate, avp{a) 1 is non- 
negative. In almost all the situations to be considered 
here, the second term on the right-hand side of Eq. (l2.5[) 
is completely negligible compared to the first for positive 
a. 

Note that this minimal model already addresses some 
of the concerns expressed by Cottrell. Our use of the 
Taylor stress in Eq. (|2.3p recognizes that the depinning 
mechanism is a collective effect, depending on the average 
spacing of the extended objects in the Cottrell-Nabarro 
"bird's nest." Moreover, our assumption that the aver- 
age velocity v(a) is accurately determined by only the 
depinning time Tp is consistent with Cottrell's "weakest 
link" picture. The speed at which a dislocation segment 
jumps from one pinning site to another is increasingly 
unimportant the greater that speed becomes; the slowest 
process - depinning in this case - is always the dominant 
one. It is this property of the system that gives us lati- 
tude in choosing which of various competing mechanisms 
determines the barrier height Tp. 

To make the depinning argument more precise, sup- 
pose that the speed of an unpinned dislocation, say 
is linearly proportional to the Peach-Koehler 

force: 



bo 

vdKP) = , 



(2.7) 



U P {a) = k B T P e- a / aT , 



(2.3) 



where r\ is a drag coefficient with the dimensions of stress. 
The motion of unpinned dislocations can make an appre- 
ciable difference in the total strain rate only if the drag 
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time td = 0* /vd is at least as long as the pinning time 
Tp, that is, if I* rj t q / b a > To/ /p(c). In principle, this 
inequality might be satisfied for stresses that are so large 
that the activation rate saturates, fp(cr) — > 1, but that 
are still less than rj. As will be seen, however, fp(a) 
remains small because, at high strain rates, the density 
of dislocations p, and thus the Taylor stress ut, become 
large. Therefore, it is not necessary to invoke a high- 
stress, viscosity dominated regime, even at the highest 
strain rates reported in PTW. On the other hand, the 
depinning argument, and the minimal model itself, do 
break down in Zaiser's fl2j limit of very low temperature 
and vanishingly small strain rates. At T = 0, there must 
be a finite stress that can overcome the pinning barrier, 
and the zero-temperature limit of Tp cannot be infinite 
as is predicted by Eq. (l2.6[) . As Zaiser has shown, disloca- 
tion dynamics in that regime is very different than that 
described here. 

For all of the analysis that follows, we assume that we 
are not in the extreme low-temperature regime, so that 
v = vp(a). We also assume that the mean free path i* 
is proportional to t, and simply let I* — t. (A numerical 
factor here can be absorbed into the definitions of other 
parameters.) Then we write 
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(a, p) = e? 1 t q = by/p\f P (a) - f P (-a) 



(2.8) 



The dimensionless plastic strain rate q(a, p) compares 
rates of mechanical deformation with intrinsic micro- 
scopic rates such as atomic vibration frequencies. Or- 
dinarily, laboratory stress-strain curves are measured at 
q <C 1. However, the highest rates in PTW approach 
q ~ 1, as do strain rates in fracture or in sheared granu- 
lar materials. We discuss the high strain rate regime in 
SeclVlIl 

Before going any further, it is useful to solve Eq. ([2.8p 
for the ratio o j ot- Except in the elastic region, at vanish- 
ingly small strains, the total strain rate is entirely plastic 
to a very good approximation, and ct/ctt is easily large 
enough that the reverse-stress term on the right-hand 
side of Eq. (|2.8p is negligible for positive stress. Dropping 
that term, we find 



In 
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u(T,p,q) (2.9) 



Note that q appears in Eq. (|2.9|) only as the argument 
of a double logarithm; thus rate hardening for q <C 1 
is extremely slow, as is observed experimentally. Most 
importantly, v(T, p, q), is a very slowly varying function 
of all its arguments, indicating that a is always a slowly 
varying multiple of the Taylor stress within the range of 
validity of this approximation. 



EFFECTIVE-TEMPERATURE 
THERMODYNAMICS 



III. 



The constitutive relation in Eq. (l2.8[) must be supple- 
mented by a theory that, ultimately, will tell us how to 



evaluate the dislocation density p that appears there both 
explicitly and as the argument of the Taylor stress &t- 
The thermodynamic ideas that we use to develop this 
theory are best - in fact, necessarily - expressed in terms 
of the effective temperature x defined in Eq. (|l.ip . These 
ideas are discussed in detail in [13] ; the results presented 
in [2(| are especially relevant to what follows. For sim- 
plicity, we focus here only on the aspects of the effective 
temperature theory that are relevant to the dislocation 
problem. 

We start by assuming that the internal degrees of free- 
dom of a solidlike material can be separated into two, 
weakly interacting, subsystems. The first, configura- 
tional, subsystem is defined by the mechanically stable 
positions of the constituent atoms, i.e. the "inherent 
structures" of Stillinger and Weber. (2ll. l22l| The second, 
kinetic- vibrational, subsystem is defined by the momenta 
and the displacements of the atoms at small distances 
away from their stable positions. 

In general, the kinetic- vibrational degrees of freedom 
are fast variables; they relax to equilibrium on atomic 
time scales. As a result, they equilibrate rapidly with a 
heat bath, and are always at the bath temperature T . On 
the other hand, the atomic rearrangements that take the 
configurational subsystem from one inherent structure to 
another are relatively slow and/or infrequent. For exam- 
ple, the generation of dislocations by Frank- Read sources 
is an extremely slow process in comparison with atomic 
vibration frequencies. The crucial point is that, when the 
dynamic coupling between the two subsystems is weak - 
as is the coupling between dislocations and thermal fluc- 
tuations - then the configurational degrees of freedom 
can be out of equilibrium with the kinetic-vibrational 
ones. In that case, under circumstances specified more 
carefully in the configurational subsystem has an 
effective temperature of its own, denoted here by \. 

For crystalline solids, the inherent structures are spec- 
ified by the populations and positions of dislocations, 
grain boundaries, and other defects. The energy Uc in- 
troduced in the paragraph preceding Eq. (jl.ip is most ac- 
curately defined as the energy of an inherent structure. 
Let Af(Uc) dUc be the number of such structures with 
energies in the interval dUc , and define the dimensionless 
entropy to be Sc = In (N/No), where No is an irrelevant 
normalization constant. Then \ = dUc/dSc has the 
dimensions of energy. 

Because x is a thermodynamic temperature in the 
usual statistical sense, the most probable density of dislo- 
cations at given x is proportional to a Boltzmann factor 
in the limit where this density is small, that is, in the 
limit in which the probability of any given lattice site 
being occupied by a dislocation is much less than unity. 
Therefore, 

1 



p ss (x) 



- en/x 



(3.1) 



where a is a length scale of the order of atomic spac- 
ings, and eo is a characteristic formation energy for dis- 
locations - the energy per unit length of a dislocation 
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multiplied by some mesoscopic length such as the cir- 
cumference of a dislocation loop or a typical grain size. 
(See @, Hol ] for a discussion of how this familiar formula 
emerges in a nonequihbrium context.) The superscript 
"ss" means that p = p ss (x) only when the system has un- 
dergone enough deformation that the dislocation density 
has reached a steady-state quasi-equilibrium at the cur- 
rent value of x- Note that, in writing Eq. (l3.1[) . we already 
are invoking the second law of effective-temperature ther- 
modynamics. 

The first law of thermodynamics for this system is 



x = W - W S + Q. 



(3.2) 



Here, the left-hand side is the rate of change of the config- 
urational heat content per unit volume V. On the right- 
hand side of Eq. (|3.2[) , W — e pl a is the rate at which 
inelastic external work is being done by the stress <r; and 
Q is the (negative) rate at which heat is flowing from 
the kinetic-vibrational subsystem (the heat bath) into 
the configurational subsystem. Ws is the rate at which 
the configurational internal energy Uc is increasing, for 
example, by formation of new dislocations, at fixed con- 
figurational entropy Sc- This term is sometimes called 
the time derivative of the "stored energy of cold work," 
but that expression can be misleading. Note that nei- 
ther x dSc nor Ws dt are exact differentials; neither the 
configurational heat content nor the stored energy are 
state functions, and both incremental quantities include 
energies associated with dislocations and other internal, 
configurational degrees of freedom. 

Equation (|3.2[) becomes an equation of motion for the 
effective temperature x when we write it in the form 



c eff X 



W - W v 



(3.3) 



where c e ^ is the dimensionless, effective specific heat at 
constant volume and constant dislocation density p. The 
quantity W x is the rate at which Uc is increasing at fixed 
X instead of, as in Ws, at fixed Sc- The thermodynamic 
distinction between Ws and W x was not made correctly 
in Q and [Ioj |. It will not be important for present pur- 
poses; but there are other situations in which it becomes 
essential. (We thank K. Kamrin for pointing out this 
earlier mistake to us.) 

To evaluate the terms in Eq. Q3.3l) . we make two key 
observations. First, because of the very large energies 
associated with creation and annihilation of dislocations, 
ordinary thermal fluctuations are completely ineffective 
in driving those processes. Thus, the only scalar rate 
factor available to us is the rate W at which work is 
being done by the driving force. This quantity is non- 
negative because, in the absence of thermal fluctuations 
large enough to anneal out the dislocations or induce 
thermally assisted strain recovery, the strain rate must 
always be in the same direction as the stress. As a result, 
we may think of W as being proportional to the non- 
negative strength of mechanically induced noise in an 



otherwise quiet environment. It follows that each of the 
terms on the right-hand side of Eq. p.3[) is proportional 
to W. 

Second, we know that, in steady-state deformation at 
dimensionless strain rate q, the effective temperature % 
must approach a stationary value, say, Xss{q)- This func- 
tion has been measured directly in molecular-dynamics 
simulations of glassy materials by Liu and coworkers. 

(See also Q.) In the limit q < 1, where the 
shear rate is much smaller than any intrinsic rate in the 
system, they found that Xss{q — > 0) = xo is nonzero, 
roughly (perhaps exactly) equal to the glass transition 
temperature. In other words, these systems reach fluc- 
tuating steady states of disorder under slow shear. The 
slower the shear, the longer a system takes in real time 
to reach steady state; but the ultimate value of Xss i n 
this limit is independent of q. Liu et al. also observed 
that Xss(q) rises with increasing q and appears to di- 
verge as q approaches unity. We will need a model for 
that behavior when we study the limit of very high strain 
rates in Sec I VIII For the present, we simply assume that 
this function exists and is nonzero throughout the range 
< q < 1. 

Return now to the right-hand side of Eq. (|3.3l) . The 
preceding argument about the rate factor implies that 
W x is proportional to W. We also know that the pro- 
portionality factor must be non-negative and less than 
or equal to unity, because only a fraction of the exter- 
nal work is converted to stored energy. Finally, we know 
that W x vanishes in steady state, where all of the work 
done is dissipated as heat, W — —Q, and where x — Xss- 
Thus, for values of x n °t too far from Xss, we can make 
a linear approximation and write 



W-y 



(1 - «) 1 



X 

Xss 



W, < K < 1, 



(3.4) 



where k is a system-specific parameter that is propor- 
tional to the fraction of the external work that is con- 
verted into configurational heat. The parameter k will 
play a central role in the discussion of strain hardening 
in SecHVl 

Using the same reasoning, we can write the heat flow 
in the form 

X 



Q 



w 



Xss{q)' 



(3.5) 



The factor x appears here because a conventional heat- 
flow is proportional to the temperature difference ksT — 
X, and we know that ksT <C X- The factor x7^ en ~ 
sures that Xss retains its meaning as the steady state 
value of x- Note that Eq. (|3.5j> predicts that the differ- 
ential heating coefficient, usually denoted by the symbol 
fidiff = —Q/W, is simply equal to x/Xss- This formula 
looks qualitatively like the data shown in (2(|; it merits 
further investigation. 

Equation (|3 . 3[) now becomes 



1 - 



X 



Xss(q) 



(3, 
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It is convenient to rewrite this equation using the total 
strain e instead of the time as the independent variable: 



de 



q(v, p) 



-eff 



'ill 



1 



X 



Xss{q) 



5o 



er . 



(3.7) 



Equation (I3.6|) is closely related to its STZ-theory ana- 
log, Eq.(5.10) in [Til] - The detailed derivation presented 
in that paper includes the effects of ordinary thermal 
fluctuations, whose absence simplifies the present result. 



IV. STRAIN HARDENING AND AN 
EQUATION OF MOTION FOR THE 
DISLOCATION DENSITY 

We now need equations of motion for the stress a and 
the dislocation density p, relating both to the effective 
temperature x as determined by Eqs. (|3.6j) or (|3.7j) . 

To obtain an equation for the stress, assume that the 
total strain rate, e, is the sum of a plastic part e pl given 
by Eq. (|2.8[) . and an elastic part a/p, where p, is the shear 
modulus. Then, converting to total strain e as the inde- 
pendent variable, we have 



da 
Ik 



P) 



Q0 = CTq, 



(4.1) 



where q(a, p) is given by Eq. (|2.8l) supplemented by 
Eq. (|2.4p for the Taylor stress as a function of p. 

In writing an equation of motion for p, we run into 
a serious problem in the conventional literature. Much 
of modern dislocation theory is based on the so-called 
"storage-recover y eq uation" first proposed by Kocks and 
collaborators (5l.Tl6Hl8|: 



dp/dt = ki y/p - k 2 p, 



(4.2) 



where k\ and k 2 are p-independent parameters. Here, p 
is meant to be the density of just the stored (i.e. pinned) 
dislocations; and ^fp is the inverse of the dislocation spac- 
ing £, which is assumed to be the mean free path for 
mobile dislocations. The first term on the right-hand 
side of Eq. (|4.2l) is said to be a storage rate, and the sec- 
ond term is the rate of annihilation or mobilization of 
stored dislocations. Equation (|4.2j) is the starting point 
for a large part of the phenomenological literature in this 
field, including @ and [HJ]. If we say that the stress 
is always equal to the Taylor stress, and replace ^fp by 
a term proportional to a, we obtain the Voce equation 
27], which implies that er(e) rises linearly from zero and 
relaxes exponentially to a steady-state flow stress in the 
limit of large, positive e. With various modifications and 
the addition of other parameters, stress-strain curves of 
this kind can be made to fit a wide range of experimental 
data. 

One problem with Eq. (|4.2[) is that the left-hand side 
changes sign when the shear rate changes direction; thus, 
this equation violates time reversal and reflection symme- 
tries. The symmetry problem can be solved superficially 



by using only positive values of the strain, effectively by 
replacing dp/de by its absolute value; but such a math- 
ematical singularity cannot appear at this place in any 
physically well posed equation of motion. Eq. (|4.2l) . or 
the equivalent Voce equation as generalized by PTW, 
provides a phenomenological curve-fitting device; but it 
can have no predictive value of its own. 

In fact, important physics is missing. This phenomeno- 
logical approach provides no way to introduce and test 
any specific physical model such as the one proposed here 
in SecHU Moreover, although Eq. (|4.2l) supposedly de- 
scribes the evolution of p, it contains no connection be- 
tween the rate at which mechanical work is being done on 
the system and the rate at which dislocations are being 
created or annihilated. Nothing happens in this system 
unless work is being done on it, so it seems obvious that 
the equation of motion for p - in analogy to the equation 
for x ~ must contain the external forcing. 

Accordingly, we propose to discard Eq. (|4.2p entirely, 
and replace it by an equation based on the second law 
of thermodynamics and energy-conservation. The second 
law requires that p relax toward its most probable value 
p ss (x)- Therefore, for p not too far from p ss , we write 
the equation of motion for p in the form 



dp 
~dc 



aq(a, p) 



ID 



qo 



i 



p 



p ss (x) 



qa = €T Q . (4.3) 



Here dimensionless energy-conversion coefficient 

analogous to k in Eqs. Q3.6p and (|3.7p . and jd is the en- 
ergy per unit length of a dislocation. Thus, the prefactor 
in Eq. ([4.3[) is the rate at which dislocations are being 
created if a fraction k p of the work done on the system 
is stored in that form. Note that this prefactor is non- 
negative in accord with the second law. The work rate 
a q is a non-negative scalar, and go changes sign when e 
changes sign. 

Formally, Eq. ([4.3[) can be interpreted as having 
emerged from a Clausius-Duhem inequality requiring 
that the rate of entropy production be non-negative. It 
is shown in 0, [13] that such inequalities, and equations 
of motion of this form, follow directly from the princi- 
ple that the statistically defined entropy must be a non- 
decreasing function of time. (A similar equation was used 
in 28] to describe the approach to quasi-equilibrium of 
a density-like variable.) By invoking this principle, we 
eliminate the need at this point in the analysis to spec- 
ify the mechanisms by which the system achieves stead 
state, i.e. "dynamic recovery." In Cottrell's terms 
we can "interpret" this process as being achieved by a 
balance between many dislocation creation, annihilation, 
and transformation mechanisms; but, for "predictive" 
purposes, we need only a few parameters that somehow 
contain all of this mechanistic information. We already 
have seen that the activation energy ks Tp in Eq. (|2.3[) 
can describe a wide variety of rate-limiting processes. In 
Eqs. (|3.6[) and (|4.3p . the energy-conversion coefficients k 
and hip circ playing similar roles. 

According to the discussion preceding Eq. (|3.6p . k is 
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proportional to the fraction of the work of deformation 
that is converted into configurational heat. We expect 
k to be a temperature dependent quantity. At low T, 
the work of deformation goes primarily into producing 
new dislocations, which means that it is stored in the 
form of recoverable energy, and, according to Eq. (|3.4[) . n 
is small. At higher temperatures, where pinning forces 
are weaker and the system can explore a larger range of 
configurations, more entropy is generated, and thus more 
of the energy is converted into configurational disorder. 
Therefore, we expect k to increase with increasing tem- 
perature. As yet, however, we have no way to compute k 
from first principles; therefore, we use it as an adjustable 
parameter. 

The conversion coefficient k p necessarily contains more 
structure, because it governs, not just the rate at which 
p approaches its steady-state equilibrium value, but also 
the initial hardening rate. The onset of hardening is gov- 
erned by Eq. (|4.1|) . which is a stiff differential equation 
because the shear modulus p is about two orders of mag- 
nitude larger than the Taylor stress. The crossover from 
the initial elastic behavior, where a = p e, to the onset 
of hardening occurs when q becomes nearly equal to go- 
This happens at a value of the strain that is of the order 
of 10~ 5 in the examples to be described in Sec lVD Be- 
yond that strain, q remains essentially constant at qo - 
the strain rate becomes entirely plastic - and the stress 
required to maintain the fixed q = qo grows because the 
density of dislocations is growing in accord with Eqs. (|3.7p 
and (14. 3|) . Thus, hardening is a slow and very nearly 
steady-state process. To a good approximation, we can 
use Eq. (|2.9|) . with v(T, p, q) replaced by ^(T, p, q ), to 
deduce that the ratio a /ax is equal to a weakly temper- 
ature and strain-rate dependent constant throughout the 
hardening process. 

Now consider the onset of hardening, where the total 
strain rate is changing from elastic to plastic, and the 
plastic deformation is just becoming visible on a stress- 
strain graph. Experimental evidence (see, for example, 
Fig.21 in [|| or Fig.3 in PU) indicates that the initial 
hardening rate, in units of the shear modulus, is roughly 
independent of both temperature and strain rate, and 
has a magnitude 



6 _ f (da 
P 



p \de 



1 

20' 



(4.4) 



The question is how to incorporate this experimentally 
observed, apparently universal, onset condition into the 
equation of motion for p. One possibility might be to 
supplement the activated depinning rate in Eqs. (|2.5[ ) and 
(I2.6|) by some temperature and strain-rate independent 
mechanism; but that procedure would move us away from 
our strategy of exploring only the minimal model intro- 
duced in Secjlll For the present, therefore, we reserve 
the possibility of alternative dynamical mechanisms for 
later investigation, and simply assume here that the on- 
set physics is contained in the conversion coefficient k p . 
To implement this assumption, look at Eq. (|4.3p in the 



case where the system is beyond onset but p is still much 
less than p ss . In this case, we know that q = q$, and that 
a = isqctt, where i>o = v{T, p, qo). Therefore 



n p v a T _b k p vq p T 



Id 



Id 



Solving for p, we find 

da ^ v% p\ b 2 
de 2j D 



(4.5) 



(4.6) 



If we set the right-hand side of this relation equal to Go 
and use the approximation in Eq. (|4.4l) . the prefactor on 
the right-hand side of Eq. (I4.3[) becomes 



JD 



2 6 



(v p T b) 2 10 {vqpt b)^ 



(4.7) 



and we have completely determined all the parameters in 
Eq.(g3D. 

This procedure evades the question of why or whether 
Go/ p can be a universal ratio, independent of strain rate 
or temperature, for a wide range of experimental situa- 
tions. In this connection, note that - if— we assume that 
onset always occurs when a is some fixed multiple of ar, 
and that a fixed fraction of the work of deformation al- 
ways is converted into new dislocations at this point, then 
energy conservation implies a strain-rate- independent re- 
lation of the form of Eq. (|4.5p : 



ax 
oc — 

onset 



PT 

ID 



Thus, 



6o 



Mr 



(4.8) 



(4.9) 



which seems to be a plausible, temperature-independent 
estimate of G>o/p if p, 7d, and pr all scale with tem- 
perature in the same ways. More generally, however, we 
expect that O is sensitive to many dynamical details and 
also to sample preparation, and therefore is not actually 
an intrinsic property of a material. 



V. PARAMETERS AND SCALING 

Our next step is to identify conveniently rescaled vari- 
ables and estimate values of parameters that emerge from 
this process. 

Because the formation energy eu in Eq. (13.I[) is large, 
and the effective temperature x is t ne only energy in the 
theory that is comparable to it, we transform to the di- 
mensionless ratio x = x/ e D, and use this variable in the 
equations of motion. There are several ways to estimate 
the scale of % by purely geometric considerations. Start 
with Eq. (|3.I[) for the density of dislocations p ss , and note 
that the length scale a is the average spacing between 
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dislocations when x ~ * 00 ■ There is nothing unrealistic 
about the concept of an infinite x; it describes a state 
of maximum disorder in which any area a 2 is as likely 
to contain a dislocation as not to contain one. (Infinite 
"spin temperatures" are commonly used to describe mag- 
netic systems in which the moments are equally likely to 
be aligned parallel or antiparallel with some axis.) Arbi- 
trarily large values of x play prominent roles in the high 
strain-rate analysis to be discussed in Sec lVIIl The den- 
sity a -2 , where x ~ ► 00 1 may be the value of p where the 
interactions between dislocations are energetically com- 
parable to their formation energies, so that the Boltz- 
mann approximation breaks down. Thus, it seems rea- 
sonable to guess that a might be about ten atomic spac- 
ings. 

Now recall that xo = Xss(q — > 0), via Eq. (|3.1[) . sets 
the density of dislocations when the system undergoes 
arbitrarily slow deformations for arbitrarily long times. 
We propose that this definition of xo be interpreted 
very roughly as a system-independent geometric crite- 
rion, weakly analogous to the idea that amorphous ma- 
terials become glassy when their densities are of the or- 
der of maximally random jammed packings, or to the 
Lindemann criterion according to which crystals melt 
when thermal vibration amplitudes are of the order of a 
tenth of the lattice spacing. In that spirit, we guess that 
Xo is the dimensionless effective temperature at which 
the spacing between dislocations I is roughly ten length 
scales a, or about one hundred atomic spacings. This 
estimate is consistent with the dislocation spacings at 
the upper end of the graph of ,/p versus stress shown 
in Fig.l of Mecking and Kocks 17|. Thus we guess that 
1/Xo ~ 2 ln(10) ~ 4 and, from here on, use xo = 0.25. In 
the future, when we consider more complex models con- 
taining multiple energy scales comparable to eo, we may 
be able to obtain a better estimate of xo- For the mo- 
ment, however, we presume that any inaccuracy in this 
estimate is compensated by variations in other parame- 
ters such as the ratio b/a. 

This rescaling of the effective temperature suggests 
that we define 



(5.1) 



Then we rewrite Eq. (|2.8p in the form 

q(a,p) = eP l r = (b/a)q(a,p), (5.2) 

where 



<?0, P) 



fp(v) - fp(-v) 



(5.3) 



Equation 



becomes 



JL«lnf^|-h. 



(XT 

The Taylor stress is 

(X T = P>T VP 



2 U 2 



= u(T,p,q). (5.4) 
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FIG. 1: (Color online) Stress-strain graphs for Cu at relatively 
high temperatures as shown. The data points are taken from 
PTW [l5[, Fig.2. Temperatures and strain rates are shown on 
the graph. The theoretical curves are computed with Tp = 
40, 800 if and X ss = Xo = 0.25. For T = 1173 if, p T = 1343 
MPa, and fC = 0.11 (MPa) -1 . The initial values of x are 
X in = 0.18 for the upper curve with large strain rate, and 
Xin = 0.22 for the bottom curve, where the small strain rate 
apparently allows Xm to be close to its steady-state value. 
For T = 1023 K, p T = 1490 MPa, K = 0.055 (MPa)" 1 , and 
Xzn = 0.185. 



We also use the definition of q in Eq. (|5.3[) to rescale the 
steady-state effective temperature: 



Xss(q) = Xss{q)/e D - 



(5.6) 



p T = (b/a)p 



T- 



(5.5) 



Using q instead of q as the dimensionless measure of 
plastic strain rate means that we are effectively rescal- 
ing To by a factor b/a. For purposes of this analysis, 
we assume that (a/6) To = 10~ 12 sea, independent of 
temperature; and we use q = 10~ 12 e pl for converting 
from q to measured strain rates. This estimate of To is 
about the same as the atomic vibration time used by 
PTW. The dimensionless rate q is approximately equal 
to unity at the upper edge of the PTW data; and ?«1 
throughout the thermal-activation region. It follows that 
Xss = Xo — 0.25 in the latter region. 

We now use the steady-state flow stresses from the 
high-temperature stress strain curves for Cu shown in 
PTW Fig.2 - shown here in FigQ]- to compute Tp, which 
we assume to be independent of T. We also use this infor- 
mation to evaluate the modulus px that appears in the 
definition of the Taylor stress in Eq. (|5.5l) . In Eq. (|5.4[) . 
we set p — p ss — exp (— 1/xo) with xo = 0.25. Using just 
the steady-state flow stresses for the two curves at T — 
1173 if, with dimensionless rate factors q — 9.6 x 10~ 10 
and 6.6 x 1(T 14 wc find T P = 40, 800 K and p T = 1343 
MPa. At the lower temperature, T = 1023 A", we find 
p T = 1490 MPa. At a yet lower temperature, T = 298 K, 
using data shown in Fig[21 we find px — 1600 MPa. 
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FIG. 2: (Color online) Stress-strain graphs for Cu at T = 
293 K and two different strain rates as shown. The data 
points are taken from [Ts| and [2S|. For both curves, pr = 
1600 MPa, K. = 0.007 (MPa)" 1 , and % n = 0.18. 



VI. STRAIN HARDENING: NUMERICAL 
EXAMPLES AND COMPARISONS WITH 
EXPERIMENT 

In the scaled variables, the equations of motion that 
we need for the strain-hardening analysis are: 



da 



1 



d-X r gfo P) 

—r~ = r^a — ; 

de q 



5o 



1 - 



Xss{q) 



q a = (a/b) er ; 



K = 



.1 



c e ff e D 



and 



K,p a q(a, p) 



dp 

de h>(T,p,q ) 2 q 



p ss (x) 



where 



26 p 



p T 10 p T 



(6.2) 



(6.3) 



(6.4) 



JC and IC P are temperature dependent quantities (the lat- 
ter via p and Pt), with the dimensions of inverse stress. 

The experimental data sets on which we base our 
strain-hardening analyses are shown in Figsfl] and [5] 
Both figures show constant strain-rate, stress-strain 
curves for Cu, the first (taken from PTW [H|, Fig.2) 
for two relatively high temperatures, the second (from 
and [1^1) at about room temperature. Note that, for 
both T = 1173 K in Fig H] and T = 298 K in FigH the 
two strain rates shown differ by factors of about 10 6 . 

The theoretical curves in both figures have been com- 
puted by integrating Eqs. (|6.1l) . (|6.2I) . and (16. 3[) . In 



FIG. 3: (Color online) Transition between elastic behavior 
and the onset of hardening at very small strains, for the five 
stress-strain curves shown in Figs[T] and [21 From bottom to 
top: T = 1173 K, e = 0.066 sec -1 (black); T = 1173 K, e = 
960 sec" 1 (green); T = 1023 K, e = 1800 sec -1 (red); T = 
298 K, e = 0.002 sec" 1 (cyan); T = 298^, e = 2000 sec" 1 
(blue). 



Eq. ([6.2p . we have set x ss (q) = xo because we consider 
only q < 1. We have set 6 = p/20 as in Eq. ([rl4l . 
which is consistent with the data shown in FigF2]and, in 
the absence of small-strain data in FigJTJ provides a plau- 
sible extrapolation to zero strain. We have assumed that 
p = 50 GPa at T = 298 K, and that pr scales like p as a 
function of temperature. Thus, at all temperatures, we 
have used the low-temperature ratio p = 31 pr and, in 
evaluating the right-hand side of Eq. (|6.4l) . have written 
K p = 3.1/pt- 

For initial conditions, at e = 0, we have used cr(0) = 0, 
and have arbitrarily chosen a non-zero initial value of the 
dislocation density, p(0) = 10" 7 , because we are not es- 
pecially interested in the earliest stages of deformation of 
ultra-pure crystals. The initial value of p determines the 
stress at which the onset of hardening occurs, but seems 
to have little or no effect on the subsequent hardening 
so long as p(0) <§C p ss . Accordingly, the only parame- 
ters that we have varied freely to fit the data (having 
determined values of pr from steady-state flow stresses) 
are IC and the initial value of x, say, Xin- The values of 
these parameters are shown in the figure captions. As 
expected, K is independent of strain rate, and decreases 
substantially with decreasing temperature. 

The elastic regions near the origin are invisible in both 
figures. To show what is happening there theoretically, 
we zoom in to very small strains in Fig[3j where we see 
the expected elastic behavior, a = pe, breaking off to 
the onset of hardening at stresses in the range 0.1 — 1.0 
MPa. 
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FIG. 4: (Color online) Dimensionless hardening rate Q(a)/p 
for the five stress-strain curves shown in Figs[T] and [2] The 
color code is the same as in FigO going left to right from 
higher to lower temperatures, and from smaller to larger 
strain rates. 

The hardening rates, 

m=i* (6 .5) 

p fi de 

are shown in Fig |4] for all five of the cases shown in FigsQ] 
and El At the lowest temperature, T — 293 K in FigH 
and the two right-most curves in FigHl the theoretical 
fits are almost insensitive to the initial effective temper- 
ature, Xi n . Here, hardening occurs in two stages roughly 
corresponding to stages II and III defined by Kocks and 
Mecking in Q. For a while after onset, the dislocation 
density grows in a deterministic way, independent of the 
effective temperature x, and is governed only by Eq. (|6.3j) 
for p <C p ss (x)- As is obvious both in the data and the 
theory, the early-stage hardening quickly becomes sen- 
sitive to strain rate and temperature; only the initial 
rate Qq is a strain-rate independent constant. A later 
stage of hardening, starting near the inflection points on 
the right-most curves in Fig|4j sets in when p becomes 
comparable in magnitude to p ss {x)- Then, hardening 
begins to be controlled by changes in x as determined 
by Eq. (|6.2l) . By this point in the process, however, Xin 
is largely irrelevant, and only the rate at which x ap- 
proaches its steady-state value xo is important. 

The situation is quite different for the high- 
temperature behavior shown in Fig(TJ and by the three 
left-hand hardening curves in Figj4] Here, p grows 
rapidly toward p ss (x)> and the dynamics of x is already 
important at the earliest data points shown. Since we 
have no data closer to onset, we are free to vary both K 
and Xi n , subject only to the constraints that the theoreti- 
cal curves join smoothly to the later-stage data, and that 
K, should be independent of strain rate. Our best-fit pa- 
rameters seem to make the theory look as uninteresting as 
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FIG. 5: (Color online) Effect of decreasing the initial value of 
X- All theoretical parameters are the same as those for the 
curve in Fig [T] for T = 1173 K, e = 960 sec." 1 , except that 
the values of Xin are, from left to right, 0.18, 0.10, and 0.07 
respectively. For reference, the original experimental data is 
indicated by the blue circles. 



possible, but a wider range of behaviors is consistent with 
the data. In fact, O(cr) is not a monotonically decreasing 
function if Xin is small enough but is compensated by a 
larger value of JC. 

The wide range of hardening behaviors already emerg- 
ing in this minimal model tells us that there is a rich 
variety of phenomena that can be predicted by analy- 
ses of this kind. It must be possible to use the methods 
presented here to write the equations of motion for any 
well posed dynamical model of dislocations in a thermo- 
dynamically consistent way. As mentioned earlier, this 
procedure will involve introducing separate density func- 
tions for different kinds of dislocations. 

However, even in our minimal model, there are many 
phenomena that can occur. For example, we show in 
Figj5] what happens to the stress-strain curve for T = 
1173 K. e = 960 sec. -1 , if we fix K but decrease Xin from 
its original value of 0.18 to 0.10 and then to 0.07. The 
resulting curves look qualitatively like those shown in 
Fig.l of [5| for Cu single crystals oriented at different 
angles relative to the shear stress. The similarity could 
be a coincidence; but it also could mean that the rele- 
vant populations of dislocations at each orientation have 
different formation energies and play different roles in 
the flow of energy and entropy through the system. The 
thermodynamic analysis in Sec lIIH especially the relation 
between heat generation and the effective temperature in 
Eq. (|3.5[) . might help to sort out these possibilities. 
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VII. VERY HIGH STRAIN RATES 

We turn finally to the regime of high strain rates, 
specifically to the points derived from strong-shock data 
shown in Figs. 7 and 9 in PTW for steady-state values of 
e in the range 10 9 to 10 12 sec -1 . With t ~ 10~ 12 sec _1 , 
our dimensionless strain rates q are no longer much 
smaller than unity, and therefore we must pay attention 
to the q-dependence of Xss(q)- This issue has emerged 
prominently in studies of glassy systems, and has been 
investigated in detail in [20j • Much of the following dis- 
cussion is based on the latter paper. 

To estimate Xss(q) for large values of q, note that the 
inverse function, say q = Q(x\ss), is a rate expressed as a 
function of a temperature. As such, it can be compared 
to more familiar properties of disordered systems such as 
temperature dependent viscosities or diffusion constants. 
For example, Liu et al 0, [24| found that their analog 
of Xss (q) increases rapidly with increasing q, and shows 
signs of diverging near q = 1. An Arrhenius plot of the 
low-temperature data in this range, shown in Fig. 1 of 
1201 , reveals that 



QiXa.) 



(7.1) 



where ba = e^/eo, and is an Arrhenius activation 
energy. For the glass simulations, is somewhat larger 
than the activation energy deduced directly from viscos- 
ity measurements (the analog of eo). The Arrhenius fit 
to Eq. (|7.ip in [20| is much cleaner than is the correspond- 
ing fit to the simulated viscosity as a function of the bath 
temperature. 

It is tempting to speculate that there exists a univer- 
sal, Arrhenius-like, steady-state relation between the ef- 
fective temperature and the strain rate, and that this 
relation is valid for a wide class of athermal disordered 
systems including both glass formers and crystals with 
high densities of defects. This speculation, of course, re- 
quires better theoretical justification and experimental 
tests before it can be accepted. Note, however, that it is 
supported qualitatively by the high-strain-rate, steady- 
state data of PTW. If we assume that the steady-state 
stress is always proportional to the Taylor stress, and use 
Ea.(|7.ip to evaluate p ss ~ exp (— l/xss) in Eq. (|5.5[) . we 
find 



P = 



2e A ' 



(7.2) 



The PTW results imply that ep /ca ~ 0.5, which is some- 
what smaller than the value of this ratio (~ 0.67) found 
for glass simulations in [20j], but is not qualitatively in- 
consistent with it. In short, the effective temperature 
theory predicts that rate hardening is rapidly acceler- 
ated by growth in the density of dislocations when the 
strain rate becomes comparable to the intrinsic inverse 
time scale t^ 1 . 

To complete the theoretical development, we need an 
expression for x ss (q) that crosses over from Xss — Xo 



in the small-g limit to the activated behavior shown in 
Eq. (|7.1j) . Note that the easy-to- understand saturation of 
Xss at xo in the small-g limit corresponds to the familiar 
- yet mysterious - divergence of the time scale at the 
glass transition. In (2pj and elsewhere, an approximation 
for the rate factor in the whole range, xo < Xss(q) < oo, 
is written in the form 



Q(Xss) ~ exp 



a(Xss) 

Xss 



(7.3) 



where a(x ss ) is a Vogel-Fulcher function that has been 
modified so that it cuts off smoothly near Xss — Xa, 
i.e. at the high-temperature end of the super- Arrhenius 



region: 



a(Xs 



Xi 



Xs 



— exp 

Xo 



Xs 



Xo 



XA - Xo 



(7.4) 



(The factor 3 is an arbitrary fitting parameter.) This 
form of the rate factor emerges from the excitation-chain 
theory of the glass transition [3(| HH . A more physically 
motivated expression for a(x SB ) has been introduced in 

The dashed line in FigjS] shows what happens when 
we set Q(xss) = q, solve Eq. (|7.3l) for Xss(q), and insert 
the result into Eq. (|5.4p with p = p ss = exp (— l/Xss) 
to compute the steady-state a as a function of q. The 
experimental points shown are taken directly from PTW 
Fig. 7, for which the nominal temperature is T = 300 A". 
We also include the single point at T = 298 A, e = 2 x 
10~ 3 sec -1 (the lower curve in FigJ^]) in order to show the 
theoretical connection across the entire range of strain 
rates. The best fit occurs at the PTW value, ep/2eA = 
0.25. However, the actual slope of the dashed line in FigJS] 
is closer to 0.3, meaning that the approximations leading 
to the simple estimate in Eq. (|7.2[) - especially dropping 
the q and p = p ss dependences on the right-hand side of 
Eq. (l5.4p - were not completely accurate. 

Note that ep/eA is our only adjustable parameter; the 
quantities that determine a(xss) in Eq. (|7.4l) are irrele- 
vant in this large-<7 region, where Q(x S s) is dominated 
by the Arrhenius term. Thus, our choice of ep/eA de- 
termines, not only the slope of the curve at high strain 
rates, but also the location and shape of the crossover 
between the two different rate-hardening regimes. There 
is very little flexibility in this fit. The PTW estimate for 
P was obtained, as shown in their Fig. 9, by drawing a 
straight line through data whose slope is not constant. 
In fact, our predicted slope of approximately 0.3 is cor- 
rect for the first part of the high strain rate regime, but 
the experiments indicate that the rate-hardening effect 
weakens as the strain rate increases. 

One possible explanation for this weakening is a ther- 
mal effect. The activation rate introduced in Eq. ([2.6[) 
is a strongly temperature dependent quantity; small in- 
creases in T increase this rate substantially. According to 
Eq. (l3.5[) . heat flows from the deforming configurational 
subsystem to the kinetic-vibrational subsystem at a rate 
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FIG. 6: (Color online) Stress versus strain rate for Cu at 
T = 300 K. The data points are taken from PTW [H, Fig.7. 
The solid theoretical curve is computed using e_o /2 ca = 0.25, 
Xi = 0.1, xa = 0.3, and a thermal conductivity K = 0.4. The 
dashed curve is computed without the thermal effect, i.e. with 
infinite K. 

— Q equal, in steady state, to the rate at which work is 
done in driving the deformation. If we suppose that the 
thermal conductivity between the configurational subsys- 
tem and an external heat bath is less than infinite, then 
the temperature T must increase by an amount propor- 
tional to Q. Specifically, we estimate that the heat flow 
is 

Q = qtr = K(T-T ), (7.5) 

where To is the ambient temperature and K is a thermal 
conductivity. 

The solid curve in Fig|S]has been computed by insert- 
ing Eq. (|7.5l) into Eq. (|2.9[) and solving for a as a func- 
tion of q with T = 300 if, f3 = e D /2e A = 0.25, and 
K = 0.4. The agreement is remarkably good, consider- 
ing the substantial uncertainties in both the experiments 
and the theory. Among the theoretical uncertainties is 
the fact that Eq. (17.5[) predicts the temperature T to be 
substantially larger than the melting point at the highest 
strain rates. Perhaps the system can remain in a super- 
heated solid state during the microscopically short time 
intervals over which the deformations take place. (The 
experiments use shock fronts to achieve such high strain 
rates. These are not really steady-state measurements.) 
Note also that the thermally modified theory still under- 
estimates the thermal softening effect at the very highest 
strain rates, perhaps indicating that the material is in- 
deed melting in that regime. 

To complete the comparisons with PTW, we show our 
version of their Fig. 11 in FigJTl Here, we supplement the 
300 K curve shown in FigJ5] with the stress-strain-rate 
curve for IIHK and include the experimental point de- 
duced from the flow stress at 960 sec -1 as seen in Fig|T] 
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FIG. 7: (Color online) Stress versus strain rate for three dif- 
ferent temperatures as shown on the graph. The theoretical 
parameters are the same as those used for Fig(5] 



We also draw a curve for T = 100 if (guessing that 
fix — 1700 MPa), to illustrate that strain-rate hardening 
is extremely slow in general, and becomes slower at lower 
temperatures. Even with the thermal effect included, 
there remains a substantial amount of temperature de- 
pendence at the highest strain rates. 



VIII. SUMMARY AND CONCLUDING 
REMARKS 

The dislocations in a deforming crystalline or polycrys- 
talline solid constitute a complex subsystem of the mate- 
rial as a whole. The configurations of this subsystem have 
a wide range of energies, and there are extensive num- 
bers of such configurations in any energy interval. Such 
a system has an entropy, and therefore it has a temper- 
ature. This "effective" temperature is enormously larger 
than the ordinary temperature of the solid. Although 
ordinary thermal fluctuations may affect the mobility of 
dislocations, only the work done by external forces is on a 
scale comparable to dislocation energies. Thus, to a first 
approximation, the subsystem of dislocations is decou- 
pled from the heat bath, and the effective temperature 
is a well defined property of the subsystem. The prin- 
cipal theme of this paper and its predecessors |9Hll1| is 
that any such driven subsystem is amenable to thermo- 
dynamic analysis. In particular, its equations of motion 
must be consistent with the first and second laws of ther- 
modynamics. 

We have based our thermodynamic equations of mo- 
tion primarily on two physical assumptions. The first of 
these is the role of the Taylor stress in determining the 
rate of thermally assisted plastic flow, i.e. in Eqs. (|2.4l) . 
(|2.5[) . and (|2.6p . The second assumption consists of a 
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set of relations between the effective temperature and 
steady-state properties of the system, in particular, the 
quasi-equilibrium density of dislocations in Eq. (|3.1j) . and 
the plastic strain rate in Eq. (|7.3l) . The former relation is 
simply a consequence of the fact that the effective tem- 
perature is the "true" temperature of the configurational 
subsystem. The latter has emerged recently in theories 
of glass dynamics; its form may be generally valid for all 
strongly disordered, solidlike materials. 

Our theory contains only a small number of adjustable 
parameters, all of which have physical interpretations 
in terms of known properties of materials, and most of 
which are directly measurable by steady-state experi- 
ments. These parameters are: an overall microscopic 
time scale to; a baseline value of the dimensionless ef- 
fective temperature xo that we estimate a priori from 
geometrical considerations; the height of the depinning 
barrier Tp; the shear modulus p; a reduced shear modu- 
lus pt that determines the Taylor stress; the initial hard- 
ening rate Oo/p that we assume to be temperature and 
strain-rate independent; a ratio of characteristic disloca- 
tion energies eo/eA] and a quantity K that is propor- 
tional to the fraction of the work of deformation that is 
converted into configurational heat. An analogous con- 
version factor, K. p , related to formation of dislocations, 
is determined unambiguously by ©o and other known 
parameters. Of these parameters, only JC is a temper- 
ature (but not strain rate) dependent quantity that we 
choose to fit the transient, strain-hardening curves. Our 
only other, similarly free, parameter is the initial value 
of the effective temperature, which necessarily depends 
on sample preparation. With these ingredients, we pre- 
dict strain hardening in agreement with experiment over 
a wide range of temperatures and strain rates. We also, 
with only the single adjustable parameter eo/eA, accu- 
rately predict both the crossover from moderate to very 
high strain rates and the power-law exponent j3 observed 
in the latter regime. 

This theory provides a framework for further inves- 
tigation in at least two important directions. First, as 
stated at the end of Sec lVIl we should recast it in terms 
of a set of densities for different kinds of dislocations, 
each density function appearing as an internal state vari- 
able obeying its own dynamical equation of motion. (See 
9] for a discussion of the role of internal variables in 
nonequilibrium thermodynamics.) The equations of mo- 
tion for the different kinds of dislocations can include 
the mechanisms by which these populations interact with 



each other, thus possibly modifying the simple depinning 
model used throughout this paper. In this extended de- 
velopment, we might also introduce other internal vari- 
ables such as average grain size or impurity concentra- 
tions in order to make contact with other work in this 
field. We reiterate that it must be possible to use the 
methods presented here to write the equations of motion 
for any well-posed dynamical model of dislocations in a 
thermodynamically consistent way. 

Second, we should introduce spatial heterogeneities. 
From the beginning of this analysis, we have considered 
only homogeneous systems; thus, we have ignored the 
formation of dislocation microstructures, shear-banding 
instabilities, fracture, and the like. The reason for a lack 
of steady-state data at lower temperatures and higher 
strain rates is that materials fail heterogeneously under 
those conditions. One goal of a theory such as that pro- 
posed here is to predict quantitative failure criteria. 

In fact, the dynamical equations that we have derived 
are already continuum approximations; they can easily 
be generalized to situations in which a, x anc > p are spa- 
tially varying fields. For amorphous materials, the anal- 
ogous effective-temperature theory already has provided 
a way to understand shear banding as a rate-weakening 
instability. The effective heat generated in a region that 
happens to be deforming faster than its neighbors soft- 
ens that region, accelerating the deformation rate and 
generating yet more heat. Unlike the ordinary temper- 
ature, which diffuses very rapidly, the effective temper- 
ature hardly diffuses at all. As a result, the effective- 
temperature theory predicts instabilities with realistic 
length and time scales. See [33[ for a comparison of 
this theory with molecular dynamics simulations of amor- 
phous shear banding [34| , and [H, [36| for applications to 
the dynamics of shear fracture in earthquake faults. 

An even wider range of dynamic heterogeneities is 
known to occur during dislocation-mediated deforma- 
tion. It should be a fairly straightforward exercise to re- 
peat the amorphous shear-banding analyses for polycrys- 
talline solids using the present dislocation theory. The re- 
sults clearly will be different because, unlike amorphous 
plasticity, simple dislocation mediated plasticity is rate 
strengthening. It will be interesting to learn what phys- 
ical ingredients must be added to the present theory to 
produce strain localization. An equally interesting ques- 
tion is whether some version of this theory can predict the 
diverse kinds of dynamic dislocation patterns discussed, 
for example, by Ananthakrishna. [l9| 
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